Microscopic description of the pygmy and giant electric dipole resonances 

in stable Ca isotopes. 
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The properties of the pygmy (PDR) and giant dipole resonance (GDR)in the stable '^'^Ca,'^*Ca and 
*^Ca isotopes have been calculated within the Extended Theory of Finite Fermi Sy stems {ETTVFS). 
This approach is based on the random phase approximation (RPA) and includes the single particle 
continuum as well as the coupling to low-lying collectives states which are considered in a consistent 
microscopic way. For *''Ca we also include pairing correlations. We obtain good agreement with 
the experimental data for the gross properties of both resonances. It is demonstrated that the 
recently measured A-dependence of the strength of the PDR below 10 MeV is well understood in 
our modehdue to the phonon coupling some of the strength in *^Ca is simply shifted beyond 10 
MeV. The predicted fragmentation of the PDR can be investigated in (e,e') and (7,7') experiments. 
Whereas the isovector dipole strength of the PDR is small in all Ca isotopes, we find in this region 
surprisingly strong isoscalar dipole states, in agreement with an (a, a'7) experiment. We conclude 
that for the detailed understanding of the structure of excited nuclei e.g. the PDR and GDR an 
approach like the present one is absolutely necessary. 
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I. INTRODUCTION 

The pygmy dipole resonance (PDR) or "soft" electric 
dipole resonances have attracted the interest of several 
research groups in the past decade P, Q • The PDR is 
defined as a part of the low-energy tail of the isovector 
dipole resonance (GDR) and dwells well below the GDR 
maximum. In phenomcnological models these resonances 
are described as vibrations of the excess neutrons (neu- 
tron skin) against an inert core with N=Z Therefore, 
one expects little strength in N=Z nuclei and increasing 
strength with increasing neutron excess. In exotic nuclei 
with extreme neutron excess these PDR should be espe- 
cially pronounced. As the dipole strength near the nu- 
cleon separation energy strongly influences the r-process 
in astrophysics, the PDR in such nuclei is of importance 
for the nuclear synthesis. 

Another important reason of the interest in the PDR 
are modern experimental techniques which allow to dis- 
tinguish between 1^,2+ and 1+ states In recent 
experiments new data have been obtained for several 
groups of isotopes (see 0, IE 01) including the stable 
Ca isotopes 0, ■ The latter results for the three nu- 
clei '^^Ca,'^^Ca and ^®Ca, which are the bridge between 
light and medium weight nuclei, are especially interest- 
ing. The experiment shows a large difference in the 
summed El strengths of the PDR below 10 MeV for ^°Ca 



and ^^Ca. Surprisingly the EWSR of the PDR does not 
rise with the neutron excess but is roughly the same for 
^'^Ca and ^^Ca. This contradicts the phenomcnological 
models as well as microscopic RPA calculations 

Clearly, such an analysis should be performed within 
microscopic theoretical approaches which simultaneously 
describes the PDR and GDR and other nuclear struc- 
ture properties. The models should not have free fitting 
parameters for each nucleus in order to obtain reliable re- 
sults which allow to understand the underlying physics. 
Such an analysis has been performed in ref. Q for the 
three Ca isotopes within the RPA based on the self- 
consistent density functional method and for heavier nu- 
clei within the relativistic self-consistent RPA or QRPA 
based on the relativistic Hartree-Bogoliubov model 
In both approaches the self-consistency ensured a reliable 
comparison between different nuclei. Goriely and Khan 
calculated within the self-consistent QRPA the El 
strength distributions for all nuclei with 8 < Z < 110 
between the proton and neutron drip lines using known 
Skyrme forces. In their calculation the low-lying El 
strength was located systematically higher by some 3 
MeV compared with the available data. 

In the past 15 years Kamerdzhiev, Speth, Tertychny 
and Tsclyaev have developed the Extended Theory of Fi- 
nite Fermi Systems (ETFFS) which uses the Green func- 
tion formalism, is based on the RPA and includes the 
single-particle continuum as well as the coupling to the 
low- lying phonons. It also considers the effect of the 
phonons in the ground state correlations. The authors 
have shown in numerous publications (see |lOl | and ref- 
erences therein), that this approach allows to calculate 
simultaneously low-lying as well as high lying nuclear 
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structure properties in quantitative agreement with the 
experimental data. The model has been applied to all 
closed shell nuclei. As far as the electric dipole states are 
concerned the complete spectrum has always been calcu- 
lated, but as the PDR are weak compared to the GDR 
they appear only as small fluctuations at the lower end 
of the GDR. 

The first calculations of the PDR within this model 
including pairing have been performed for Ca and Sn 
[Tl| isotopes. Here the characteristics of the PDR are 
in much better agreement with the data than the above 
mentioned self consistent calculations. This is due to 
coupling to the phonons, first of all, and due to another 
single-particle scheme. There exists a basic difference be- 
tween the present approach and self-consistent (Q)RPA 
calculations e.g. 0,013 as we will discuss below. 

Long time ago the conventional (Ip-lh) RPA and 
QRPA have been extended to include the effects of 
phonons. Soloviev ct al. introduced 30 years ago the 
Quasiparticle Phonon Model (OPM)]VA\\^. the present 
model {ET F F S) was first applied nearly 15 years ago 
in ref.0, ^3 ^"^^ around the same time Broglia et al. 
developed the Phonon Coupling Model(PCM) 0|. How- 
ever, the extended models 113 have been used only 
recently to calculate the properties of PDR in non-magic 
nuclei because an additional work had to be done to in- 
clude pairing to these models . In addition to the above 
mentioned results, the PDR calculations have been per- 
formed for ^°*P6 and for a long Sn isotopes chain 
within the QPM, for ^'^^Sn^'^Sn , ^ospfo (together with 
the giant El resonances) [l^ and for i^O,2"0,220 ^ 
within the PCM generalized to include pairing. Very re- 
cently, the El photo absorption cross sections have been 
calculated in ^^^Sn, '^'^^Sn and ^'^'^Sn with a version of 
the ETFFS (see below) |2ll| .which also considers pairing. 
So far the single-particle continuum for non-magic nuclei 
has been onl y ta ken into account in the present approach 
(ETFFS) |i,"lll|23. 

In the recent ETFFS calculations for the PDR in Ca 
isotopes @, several approximations (for details see prj l 
have been made which were responsible that the theo- 
retical results agreed only qualitatively with the data. 
The aim of the present investigation is threefold: firstly, 
we repeat the previous calculations without the approx- 
imations and apply for ^^Ca the generalized version of 
the ETFFS which includes the pairing effects, secondly 
we clarify the role of the continuum for PDR and GDR, 
which has attracted a great interest now ^3 1 ^i^d finally 
we investigate in detail the microscopic structure of the 
PDR, especially the role of the isoscalar dipole compo- 
nents. The result of the present calculation agrees well 
with the experimental data 0,1^ and it explains in a nat- 
ural way the somewhat surprising result that the strength 
of the PDR does (seemingly) not scale with the neutron 
excess. 



II. METHOD 

We have used the generalized version of the ETFFS 
with pairing which differs from the one used in [Tlj 
by taking into account the self-energy and the induced 
ph- and pp- interaction graphs on the same footing, i.e. 
by a consistent summations of the so-called g^ diagrams 
where g is the amplitude for low-lying phonon creation. 
The formalism developed for nuclei with pairing, is called 
Quasiparticle Time Blocking Approximation (QTBA), 
and is described in detail in (23i] and i21|. (However, the 
ground state correlations induced by the phonon cou- 
pling have not been included to the present calculations 
). The single-particle continuum has been included on 
the RPA level where it is taken into account correctly 
within our Green function technique in the coordinate 
representation. For details see [ifl- Therefore in our 
approach we consider the three mechanisms which cre- 
ate the width of giant resonance, namely (I) the Landau 
damping ((Q)RPA configurations), (II) the escape width 
(the single-particle continuum) and (HI) the spreading 
width (phonon coupling, or complex configurations). As 
in all our previous calculations within the ETFFS, we in- 
clude the most collective low-lying phonons, namely 3j~, 
5j" phonons for ^°Ca, and 2^,3j",5j' phonons for "^^Ca 
and ^^Ca. The collective phonons have been microscop- 
ically calculated within the RPA for ''"Ca and ^^Ca. In 
the case of ''^Ca we used the corresponding QRPA with 
pairing. The standard BCS equation have been applied 
using the particle-particle Landau-Migdal forces [2J|. In 
the present and in our previous calculations we used the 
special " forced consistency" procedure to obtain the 
energy of the spurious El state to be exactly equal to 
zero without the procedure of fitting force parameters. 

As in our previous calculations we used the effective 
particle-hole Landau-Migdal interaction: 

F(r,r')-Co[/(r) + /'(r)Ti.T2 + 

(g + .g'ri-T2)rTi-rT2]^(r-r'), (2.1) 

with the conventional interpolation formula, for example, 
for the parameter f 

fir) = fe. + ihn - fe.)poir)/po{0) (2.2) 

and similarly for the other r-dependent parameters. Here 
Pa{r) is the density distribution of the ground state of the 
nucleus under consideration and /;„ and fex are the force 
parameters inside and outside of the nucleus. The stan- 
dard values of the parameters, which have been already 
used for all the nuclei under consideration |^, 0, j are 
as follows 

/,„ = -0.002, /e, = -1.4, = 2.30, /;„ = 0.76, 

g = 0.05, g' = 0.96, Co = 300 McVfm3.(2.3) 

For the nuclear density po (r) in the interpolation formula 
we chose the theoretical ground state density distribution 
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of the corresponding nucleus, 



(2.4) 



which is more consistent than the previously used Woods- 
Saxon distribution. For that reason we had to readjust 
fex and obtained the value of fex = -1-4 used [2^. Here 
Ri{r) are the single-particle radial wave functions of the 
single-particle model used. For other details of the cal- 
culations, see [6l lTollTl| . 



III. PDR RESULTS 
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FIG. 1: The energy- weighted electric dipole strength inte- 
grated up to 10 MeV is shown for the isotopes *''Ca, ''^Ca, 
and **Ca. Dashed: experimental data taken from dash- 
dotted: quasiparticle random phase approximation, full: ex- 
tended theory of finite Fermi systems in the quasiparticle time 
blocking approximation (including pairing for '^''Ca ) 



The main results of our calculations for the gross prop- 
erties PDR are given in Fig.l and Table |I| Within the 
experimental errors, our theoretical results agree with the 
data as shown in Tablel. This holds for the total B(E1) 
strength, EWSR and the EWSR share in the energy in- 
terval up to 10 MeV. As for the PDR mean energies, 
which were defined as E ^ J:E,B^{E1)/EB,{E1), the 
agreement is rather good for '^^Ca and ^^Ca but the the- 
oretical result 8.68 MeV is too high in "'"Ca. In all the 
cases the results of the extended theory, i.e. inclusion of 
phonon coupling and single-particle continuum, give rise 
to a considerable improvement compared to the (Q)RPA 
results, and (especially for ^"'Ca) our results of ref. 0. 
The latter is due to an improvement of the theory but 
even more importantly to more careful calculations of the 
single-particle levels and of the (Q)RPA phonons, includ- 
ing the use of eq.(2.4) for them. 



In Fig.l we compare the theoretical EWSR of 
^'^Ca,^'^Ca and ^^Ca for the isovector dipole strength up 
to 10 MeV with the experimental data. The dash-dotted 
line is the (Q)RPA result which rises with neutron ex- 
cess in agreement with the phenomenological model. The 
experiment shows only a strong increase for ^^Ca com- 
pared to ^'^Ca but the strength of ^^Ca and ^"^Ca does 
not change within the experimental errors. The ETFFS 
result reproduces this somewhat surprising result. 



.a 
E 



■■= 10 



o 
u 
c 
o 

a 
o 

ID 

n 
o 
o 

£ 

a. 



RPA 
- ETFFS(QTBA) 



'"Ca PDR 



A. 



RPA 
- ETFFS(QTBA) 



Xa PDR 



7 8 9 

Ex ( MeV ) 



10 



u 

11 



FIG. 2: Pygmy resonances in ""Ca and **Ca up to 11 MeV 

The explanation is given in Fig. 2 where we compare the 
theoretical PDR-spectrum of ^^Ca and ^*Ca. Whereas 
the (Q)RPA result is in both cases around 10 MeV and 
below, the ETFFS result differs in the two nuclei. Due 
to the phonon splitting an appreciable fraction of the 
strength in *^Ca is shifted to about 11 MeV and has 
therefore not be measured in the cited experiment. If one 
adds this part ( the 10 MeV limit is somewhat arbitrary 
) the total PDR strength increase also in ETFFS with 
the neutron excess. One also realizes that the PDR for 
the open shell nucleus *'*Ca is somewhat lower compared 
to the PDR in the closed shell nucleus "^^Ca. 

In Fig. 3 the transition densities for protons and neu- 
trons of the electric dipole isovector strength up to 10 
MeV are shown for the three Ca isotopes. In ^^Ca the 
density is of purely isoscalar nature. With increasing 
neutron excess the neutron density is slightly shifted to 
a larger radius, indicating some an isovector admixture, 
but the isoscalar character remains. This is easily un- 
derstood in the framework of microscopic models. The 
isovector ph-interaction is repulsive and shifts the low- 
lying isovector strength into the region of the GDR. On 
the other hand, the isoscalar ph-interaction is attractive 
and shifts some of the high lying isoscalar dipole strength 
to lower energies and gives rise to some isoscalar collec- 
tivity below 10 MeV in contrast to the isovector case. 

In Fig. 4 the isoscalar strength distributions (with the 
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TABLE I: Properties of the pygmy dipole resonances in Ca isotopes in the (5 - 10) MeV interval. 
(The experimental data taken from [3) 







«Ca 


*4Ca 


«Ca 




(Q)RPA 


2.51 


45.6 


97.4 


lO^e^fm^ 


ETFFS(QTBA) 


3.67 


81.0 


66.2 




Exp. 


5.1 ± 0.8 


92.2 ± 15.6 


68.7 ± 7.5 




(Q)RPA 


21.48 


374.9 


897.3 


keVe^fm^ 


ETFFS(QTBA) 


31.86 


611.0 


527.7 




Exp. 


35 ± 5 


629 ± 107 


570 ± 62 




(Q)RPA 


0.014 


0.23 


0.56 


% EWSR 


ETFFS(QTBA) 


0.021 


0.38 


0.30 




Exp. 


0.020 ± 0.003 


0.39 ± 0.07 


0.33 ± 0.04 




(Q)RPA 


8.84 


8.22 


9.30 


E, McV 


ETFFS(QTBA) 


8.68 


7.54 


7.97 




Exp. 


6.80 


7.10 


8.40 




FIG. 3: Transition densities of the electric dipole strength in FIG. 4: Isoscalar electric dipole strength in Ca*° and Ca'** 
the (5-10)MeV interval 



IV. THE GIANT RESONANCE RESULTS 



operator r'^Yio) in '"'Ca and ^^Ca calculated within the 
ETFFS and RPA are shown. In ^^Ca the influence of 
the phonons are small compared to '^^Ca especially for 
the low energy spectrum. Whereas in ^"Ca the phonons 
only give rise to an energy shift one observes in '^^Ca 
also some fragmentation of the low lying strength. This 
fragmentation is due to the low-lying 2+ state which can 
be coupled to Iplh 1~ states to produce the necessary 
1~ complex configurations. The corresponding states do 
not exist in '^^Ca. The two peaks in ^°Ca at 5.8 MeV 
and 7.8 MeV are in an agreement with the experimental 
finding p7j . 



For consistency we also calculated the giant isovector 
resonances in the three calcium isotopes. Results are 
presented in Table |n] and in Figs. 5 and 6. In Table 2 
the results obtained in three different approximations are 
shown: (Q)RPA calculations with the single-particle con- 
tinuum as well as the ETFFS(QTBA) calculations with 
and without the continuum. The theoretical values are 
extracted from the calculated strength distribution under 
the condition that the three lowest energy-weighted mo- 
ments of the Lorentzian and the theoretical curve should 
coincide in the (10-40) MeV interval [H Hi]. The av- 
eraging parameter A which was used in addition to the 
single-particle continuum the giant resonance calculation 
was 400 keV. The differences between the (Q)RPA and 
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FIG. 5: Giant electric dipole resonance in Ca calculated in 
three different approximations:(I)ETFFS(QTBA) (full line), 
(II) QRPA (dash-dotted); (III) ETFFS(QTBA) without con- 
tinuum (dash-dotted) 



TABLE II: Properties of the isovector giant dipole resonances 
in Ca isotopes 





«Ca 


**Ca 


«Ca 


% EWSR 


(Q)RPA 


102.3 


106.6 


91.2 


ETFFS(QTBA)-c* 


100.4 


108.6 


91.4 


ETFFS(QTA)B 


103.7 


110.5 


93.0 


Exp. 


106.3 




119.5 


E, MeV 


(Q)RPA 


21.0 


20.4 


16.8 


ETFFS(QTBA)-c 


22.7 


23.1 


19.9 


ETFFS(QTBA) 


20.9 


21.2 


18.4 


Exp. 


20.0 




19.6 


crmax, mb 


(Q)RPA 


137.2 


146.6 


170.6 


ETFFS(QTBA)-c 


137.2 


129.2 


113.7 


ETFFS(QTBA) 


75.6 


80.9 


85.4 


Exp. 


95.0 




102.7 


F, MeV 


(Q)RPA 


2.9 


3.2 


2.6 


ETFFS(QTBA)-c 


3.2 


3.8 


4.1 


ETFFS(QTBA) 


6.2 


6.4 


5.7 


Exp. 


5.0 




7.1 



* -c stands for calculations without continuum. 



full ETFFS(QTBA) arc most striking for the width F 
and the maximum of the photo absorption cross section 
o'max, where the results differ nearly by a factor of two. 
We obtain good agreement with the experiments only 
within our full model, i.e. including the phonons and the 
single-particle continuum. Here we want to stress that 
we calculate GDR and the low-lying states simultane- 
ously within the same configuration space and with the 
same forces. 



FIG. 6: Transition densities for Ca; the isovector electric 
dipole strength is integrated in the (15-30) MeV interval 



As compared with our earlier results for °Ca and ^^Ca 
^liUS 1 here the ground state correlations induced by the 
phonon coupling were not taken into account and we have 
chosen a different energy interval for the summation. In 
addition, a different definition for the mean energy and 
another " refining" procedure have been used. Like in 
the previous calculations of giant multipole resonances in 
closed shell nuclei ^3 the inclusion of phonon coupling 
, i.e. the calculation within the full ETFFS (QTBA), 
increases the widths by a factor of 2 and (except for '^''Ca) 
shifts the mean energies by 1.0-1.5 MeV to higher energy. 

In Fig. 5 the results for the GDR in '^'^Ca are shown, cal- 
culated again within the three approximations. In the full 
model an appreciable fraction of the strength is shifted 
to much higher energies compared to the QRPA result 
which give rise to an asymmetric shape. Finally in Fig. 6 
the transition density for the GDR in "^^Ca is plotted. It 
has the expected isovector shape. 

In order to study the role of the single-particle contin- 
uum we have calculated the giant resonance (see Table 2 
and Fig. 5) and FDR characteristics without the contin- 
uum. It should be expected that the role of the single- 
particle continuum is important in light nuclei. It is of 
special interest to consider the question within the same 
calculational scheme for the Ca isotopes, including the 
non-magic '*^Ca, because there is no information about 
the role of the continuum in the approaches of similar 
kind in non-magic nuclei. One can see from Table 2 and 
Fig. 5 a very noticeable influence of the continuum for the 
integral characteristics in all the Ca isotopes as well as 
for the strength distribution. Moreover, the continuum's 
influence is noticeable even for the FDR : our full calcu- 
lations have given 0.38% of the EWSR in the (5.0-10.0) 
MeV interval for ^^Ca whereas the calculation without 
continuum gave 0.19%. The exclusion of the continuum 
also results in a considerable redistribution of the FDR 
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strength. These features can be studied in experiments 
with electrons and gamma-rays. 

V. SUMMARY 

We have investigated the low-lying (PDR) and high- 
lying (GDR) electric dipole states in the magic '^''Ca and 
"^^Ca and non-magic 44Ca within the ETFFS(QTBA) 
which, in addition to the standard (Q)RPA, takes into 
account the single-particle continuum and phonon cou- 
pling. In the present approach we corrected for the dou- 
ble counting due to the phonons in our generalized prop- 
agator. Good agreement with experiment available for 
integral characteristics has been found. 

The agreement for the PDR is much better as com- 
pared with the previous ETFFS calculations 0. The 
PDR, i.e. the isovector dipole strength below 10 MeV, 
is less than 1 percent of the energy weighted sum rule 
because the strongly repulsive isovector force shift the 
strength into GDR region. In order to investigate the 
structure of the PDR further we calculated the corre- 
sponding transition densities. It turned out that these 
densities are isoscalar with some isovector admixture in 
''^Ca and ''^Ca which is due to the strongly attractive 
isoscalar interaction which shifts the high-lying isoscalar 
strength below 10 MeV. 



In Ca we obtain two prominent isoscalar dipole states 
(induced by the external field r'^Fio) at 5.8 MeV and 7.8 
MeV. The predicted isoscalar strength can be detected 
in (a,a'7) experiments, see p^ . 

We have shown that the role of the single-particle con- 
tinuum is very noticeable for the GDR and quantatively 
important for the PDR in the Ca isotopes considered. 
The same should be true for lighter nuclei too. Thus, 
taking into account the new effects as compared with the 
standard RPA or QRPA , i.e. the phonon coupling and 
single-particle continuum , is necessary to explain the 
properties of both these resonances. 
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